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We investigate the Casimir interaction between two parallel metallic cylinders and between a 
metallic cylinder and plate. The material properties ol the metallic objects are implemented by the 
plasma, Drude and perfect metal model dielectric functions. We calculate the Casimir interaction 
numerically at all separation distances and analytically at large separations. The large-distance 
asymptotic interaction between one plasma cylinder parallel to another plasma cylinder or plate does 
not depend on the material properties, but for a Drude cylinder it depends on the dc conductivity 
a. At intermediate separations, for plasma cylinders the asymptotic interaction depends on the 
plasma wave length A p while for Drude cylinders the Casimir interaction can become independent 
of the material properties. We confirm the analytical results by the numerics and show that at short 
separations, the numerical results approach the proximity force approximation. 



I. INTRODUCTION 

Effective interactions between cylinders are an impor- 
tant parameter in synthesizing and analyzing nanomet- 
ric systems. This is due to the fact that many important 
nanostructures such as carbon nanotubes, nanowires and 
even the tobacco mosaic viruses have cylinderical shapes. 

From the perspective of the experimental Casimir force 
studies, nano-cylindrical shapes are an optimal candidate 
for precision Casimir force measurements, in comparison 
to spheres for two reasons: (i) their effective area of in- 
teraction is larger [J [5] , and (ii) mechanical oscillation 
modes of quasi-one-dimensional structures can be probed 
with high precision [3]. 

Under many circumstances, van der Waals or Casimir 
forces have the dominant contribution in the effective in- 
teractions of nanostructures, which lead to various in- 
teresting phenomena in nanosystesm. For example in 
nanomechanical devices, Casimir interaction causes stic- 
tion [4j |5] , and thus a good understanding of these forces 
leads to improvements in the design and efficiency of such 
nanosystems. In another example, Casimir interactions 
between single walled carbon nanotubes (SWCNT) with 
different chirality become important in separating a poly- 
disperse solution of SWCNT in fractions of equal chirality 

> M- 

The applicability of the Casimir interaction is not lim- 
ited to synthetic cylindrical objects. There are numerous 
examples of long macromolecular structures with cylin- 
drical shape in nature such as the tobacco mosaic viruses, 
microtubules of nagella and A-band lattice of myosin fil- 
aments in cross strained muscles [S] , and hence knowl- 
edge of the interaction between cylindrical shapes is also 
important for the biological sciences. It should be noted 
that in some biological systems composed of cylindrical 
particles which are packed in an array, the separation be- 
tween the particles can be several times larger than the 
diameter of the cylinder [7] . 



The Casimir interaction per unit length for two parallel 
perfectly conducting cylinders or a plate and cylinder at 
a separation distance d is £ / L ~ Hc/d 2 , up to a logarith- 
mic factor [HI [TO] • It decays only slowly compared to the 
retarded interaction £ / L ~ hcR 4 /d 6 between two insu- 
lating cylinders that do not support large-scale collective 
fluctuations [TT] , 

It has been demonstrated that Casimir interactions 
strongly depend on the combined effects of shape and 
material properties, see, e.g., [T!2HTT] . The interplay is 
particularly strong for quasi one-dimensional conduct- 
ing materials due to strongly anisotropic collective charge 
fluctuations. In addition, approximations of the Casimir 
force between cylinders and plates [5] have also shown 
that the temperature dependence varies based on the de- 
scription of the material properties. Thus there is a need 
for exact calculations of the Casimir force for cylindri- 
cal shapes taking into account the realistic material re- 
sponse. 

Most studies of interactions between one-dimensional 
systems over a wide range of separations concentrate 
on perfect conductors and insulators. However, low di- 
mensionality in combination with finite conductivity and 
plasmon excitations should give rise to interesting new 
effects that might be probed experimentally using, e.g., 
the coupling to mechanical oscillation modes. The often 
employed technique for these effects, the proximity force 
approximation (PFA) cannot capture the correlations of 
shape and material response since it is based on the in- 
teraction between planar surfaces. A number of stud- 
ies have been performed for the short separation regime 
mainly focused on the corrections to the Proximity Force 
Approximation [TSl - BTj . 

Van der Waals interaction between cylinders (and 
plates) have been studied for certain frequency depen- 
dent permittivities [SJ [TOJ I^IEB] . In one of the earliest 
study, the van der Waals interaction has been calculated 
between two parallel thin filaments described by one di- 
mensional (ID) plasmon and electromagnetic excitations 



[52]. This work predicts asymptotic forms of the inter- 
action energies at large separations accurately but the 
range of validity of the asymptotics remain unclear. 

In another work |23j . the Casimir interaction is ob- 
tained for conducting cylinders described by delocalized 
coupled ID plasmons at zero temperature. The response 
function of the plasmons are given by the random phase 
approximation (RPA) . The specific choice of RPA has the 
advantage that locality, additivity and R~ 6 contributions 
are not involved in the calculations. The energy is ob- 
tained by using the mode summation method, which is 
equal to the sum of the separation dependent zero-point 
plasmon modes. The Casimir energy is attractive and 
decays as ~ dr 2 apart from a logarithmic part. This 
result was later confirmed by a quantum Monte Carlo 
simulation |25j . Using the same material description and 
calculation technique, the large-distance Casimir energy 
was obtained for crossed wires with a small crossing an- 
gle. For conducting and semiconductor wires, apart from 
the logarithmic and angular parts, the large-distance in- 
teraction energy decays as ~ dr l and ~ d~ 4 , respectively 

In a completely different approach, employed for per- 
fectly conducting cylinders and plate, the Casimir energy 
is calculated for all separations using a path integral rep- 
resentation for the effective action which yields a trace 
formula for the density of states [HI [TUj . Furthermore, 
the Casimir interaction between a SWCNT and a plate 
is studied for large and short separation regimes using 
the Lifshitz formula [18]. 

The full interplay between shape and material effects 
is not transparent in the previous studies as they are 
limited either to perfect metals or to asymptotic lim- 
its. Here, we employ the scattering approach to inves- 
tigate the Casimir interaction between parallel metallic 
(circular) cylinders, and a metallic cylinder and a metal- 
lic plate. The material properties of the objects are de- 
scribed either by the plasma or the Drude dielectric func- 
tion. Some of the results have been reported in Ref. [T7] , 

The outline of this work is as follows: in Sec. [IT] we 
summarize the scattering method and the assumed ma- 
terial properties of the cylinders. 



general expression 



III we obtain 



In Sec. 

analytical results for the interaction at distances much 
larger than the cylinder radii. In Sec. [IV] the Casimir in- 
teraction is calculated numerically for different material 
properties over a wide range of separations. Section \V\ is 
the summary. 



II. METHOD 

We consider the two following systems (i) two infinitely 
long parallel cylinders, and (ii) an infinitely long cylinder 
parallel to an infinite plate. Assuming placed in vacuum, 
we calculate the Casimir interaction in these two systems 
employing the scattering formalism [11] , The Casimir 
energy of two objects at zero temperature is given by the 
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(1) 



where K is the Wick-rotated frequency and the matrix 
N factorizes into the scattering amplitudes (T-matrices) 
and translation matrices that describe the coupling be- 
tween the multipoles on distinct objects. While the mate- 
rial properties and shapes of the objects are contained in 
the T-matrices, the distance between objects is encoded 
in translation matrices. 

To implement the material properties, we consider 
plasma, Drude and perfect metal cylinders with mag- 
netic permeability \x = 1. The Drude dielectric function 
on the imaginary frequency axis is 



e(icK,) = 1 + 



(27T) 2 

(A pK ) 2 + A ct k/2 



(2) 



with conductivity a and X a = 2irc/<r. Equation (pi) re- 
produces the plasma model for A CT — > 0. 

Since the matrix N differs for parallel cylinders and 
cylinder-plate systems, in the following we describe N for 
both setups. 



A. Two Parallel Cylinders 

Consider two infinitely long, parallel cylinders with 
equal radii R and with their axes separated by a dis- 
tance d and aligned along the z-axis. The matrix N is 
diagonal in the z-component k z of the wave vector due to 
translational symmetry. The matrix elements for electric 
(E) and magnetic (M) polarizations (a, j3 — E, M) and 
partial waves m and m' are 



H 



a/3 \ A rporf 

k z mm' / j k z m 



E 



14 



12 T 7/9 7/21 
ran k z n k z nm' 



(3) 



with T the cylinder T-matrix, see Appendix |A 1| The 
translation matrix U 12 relates regular cylindrical vector 
waves to outgoing cylindrical vector waves, see Appendix 
[B] The translation matrices do not couple different polar- 
izations and for both E and M-polarization, their matrix 
elements are given by 



U? nn , = (-l) n K n _ n ,(pd). 



{-i) n - n 'ui 2 



21 



K 



(4) 



with p = yjn 2 + k 2 and K n (x) the modified Bessel func- 
tion of the second kind. 

Since N is diagonal in k z the determinant in Eq. (nl) 
factorizes into determinants at fixed k z , and the sum over 
kz moves in front of the logarithm. After taking the 



continuum limit, J^k 
length L becomes 
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|j J-oo dk z , the energy per unit 



dfc z lndet(l-N). (5) 



Here the determinant is only over the discrete partial 
wave index n. 



B. Cylinder - Plate 

Next we consider a cylinder with radius R parallel to a 
plate. We assume that the cylinder is aligned along the 
z axis and the plate is in the y — z plane. The distance 
from the center of the cylinder to the plate is d. The 
matrix N for this geometry is 



K' 



a/3 



E 



T" 7 Mf , 

k z m k z mm' 



(6) 



with 



k z mm' 

< E 

Y=E,M 



(I f< ,, 



-2d^/k^+K 2 



2v/kT 



D kzml M±l' T U ± D k ± Y,k z m'l3 (1 ~ 2( V,/3) 



(7) 



where k y is the y component of the wave vector, kj^ = 
(k y ,k z ), the matrix D n j c ^ a ] i± p converts vector plane 
wave functions and cylindrical vector wave functions, see 
Appendix C and X£ i s the dielectric plane T-matrix 
presented m Appendix |A2| The energy of this system 
can also be obtained by Eq. ^ as N is diagonal in k z , 
Note that the determinant is not related to k z ; thus, we 
suppress all the k z indices in what follows. 



The second condition ensures that the Drude dielectric 
function is large compared to one, i.e., metallic behavior 
is pronounced. At small frequencies k but fixed k z /n, for 
Drude cylinders T EE ~ k, while for plasma and perfect 



metal cylinders Tg 



EE 



1. 



Since T MM 



K , 1 Q 



EM 



T^ E = and higher order elements associated with n ^ 
scale as k 2 '™', we consider only the T EE elements at 
large separations. 

The Casimir interaction between conducting cylinders 
is intricate and no simple analytical expression that ap- 
plies to all distances can be obtained. However, using 
Eqs. ^, (3} and (J6]| along with T^ E given in Eq. (|8j), 
the asymptotic interaction at large separations, d ~^> R, 
can be evaluated in various limiting cases. 

To derive the large-distance asymptotic Casimir po- 
tential energy, we employ the identity log det = Tr log 
and expand the integrand in Eq. |5]) in powers of N, cor- 
responding to a multiple scattering expansion. The one- 
scattering approximation, sufficient for large distances, 
yields 
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(9) 



The element Noo of the N-matrix yields the dominant 
contribution to the Casimir energy at large distances 
since higher order elements involve higher powers of k. 
To this end, we will consider only this term for the rest 
of this section. 



A. Parallel cylinders 



III. LARGE-DISTANCE ASYMPTOTIC 
CASIMIR ENERGIES 

To find the asymptotic form of the Casimir interaction 
at large separations d >> R, one needs to obtain the T- 
matrix expressions for a cylinder and a plate. Using the 
dielectric function given in Eq. (pi), the asymptotic form 
of the cylinder T-matrix elements for E polarization and 
11 = at small frequencies (k <C 1, k z j k fixed) reads (see 
Eg. (|A14t), 



T, 



EE 



2r?2 



p 2 R 



C(K)-p 2 R 2 ln(pR/2) 



(8) 



where C(n) depends on the dielectric properties of the 
cylinder. For a perfect metal cylinder C(k) = and for 
a plasma cylinder C(k) « A p 2 k 2 /(27r 2 ) if the plasmon 
oscillations cannot build up transverse to the cylinder 
axis as the diameter is too small, i.e., fl « A p . In the 
opposite limit R 3> X p , we reproduce the perfect metal 
form of the T-matrix, i.e. C(k) sw 0. For the Drude 
model, C(k) = A <t k/(4tt 2 ) if k < A ff /A 2 , 1/A CT . The 
first of the two conditions implies that Drude behavior 
dominates over plasma behavior, i.e., the second term in 
the denominator of Eq. (l2| is larger than the first term. 



For two parallel cylinders, considering the fact that at 
large separations T EE yields the dominant contribution, 
the trace of the matrix N is approximated by 



Tr[N^] = T ( 



EE 7 yl2 rj^EE 7 j21 



"00 1 o 



U 



00 



(10) 



Using Eq. (10 1 in Eq. ^ and changing integration to 
polar coordinates k, = pcos(9)/d and k z = psm(6)/d, we 
obtain 



hcL 



2ir 2 d 2 ln 2 {2d/R) 

x Jo dp Jo de ^+c 1 ( P ,e))(i + C2(p > e)) 



pKlip) 



(ii) 



where C; (p, 8) describes the material properties of the 
cylinder i. For a perfect metal cylinder Ci(p,8) = 0, for 
a plasma cylinder d(p, 9)=t; cos 2 (9) with 



(12) 



2ir 2 R 2 ln(2d/R) 



in the limit R <C A p , and for a Drude cylinder Ci{p, 0) = 
£' cos(9)/p with 
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\„d 



4n 2 R 2 \n(2d/R) ' 



(13) 



For two perfect metal cylinders the integral in Eq. (11) 
can easily be calculated and yields 
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hcL 
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8nd 2 \n(2d/R) 2 ' 



(14) 



which agrees with the results in Refs. [HI HO] ■ 

For plasma cylinders with plasma wave length A p , the 
integrations in Eq. (11) yield 
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with 



/(*) = 



x + 2 

{X+ l) 3 / 2 



(15) 



(16) 



It is impo rtan t to consider Eq. ([15]) in two limiting cases 
for £, Eq. Q. In the limit £ < 1 or \n(2d/R) > A 2 /i? 2 , 
the energy simplifies to the perfect metal energy given 
in Eq. (14), i.e., the interaction between conducting 



cylinders is universal at large distances in this regime. 
However, in the opposite limit ( > 1 or equivalently 
\n(2d/R) <C X 2 /R 2 , the Casimir energy becomes 
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8y / 2A p d 2 ln 3/2 (2d/i?) 



(17) 



This shows that the universal form is applicable only 
beyond an exponentially large crossover length d ~ 
Rexp(\ 2 /R 2 )^^> R. Below this scale, and infact in any 
practical situations, the interaction is material depen- 
dent, see Fig. Qa. 

For a plasma cylinder with plasma wave length A p 
parallel to a perfect metal cylinder the integrations in 
Eq. Jill yield 
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(18) 



Similar to parallel plasma cylinders, we consider two lim- 
iting cases for £. In the limit £ <C 1 or \n(2d/R) ^> 
X 2 /R 2 , we obtain the perfect metal energy given in 
Eq. ( 14), and the conducting cylinders' interaction is uni- 



versal at large distances. In the opposite limit ( > 1 or 
equivalently \n(2d/R) <C X 2 /R 2 , the Casimir energy be- 
comes 
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(19) 



For Drude cylinders with the characteristic length A CT , 
at large separations d 3> R, A CT , we find a rather dis- 
tinct behavior that deviates from naive expectations for 
universality. In this case the integrations in Eq. (11) 



cannot be performed analytically. Therefore we first cal- 
culate the angular integral which gives a complicated ra- 
dial function, and then we expand the resulting radial 
integral for small and large £', Eq. (13). The radial in- 
tegrals can be calculated easily in these two limits. For 
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FIG. 1. Summary of the different forms of interaction be- 
tween two cylinders and a cylinder and a plate. Shown are 
the rescaled interaction energies per cylinder length, £/(hcL). 

(a) Interaction involving a plasma (P) cylinder with another 
plasma cylinder, a perfect metal (PM) cylinder or a plate. 
The asymptotic results apply sufficiently far away from the 
separating curve ln(d/_R) ~ (\ P /R) 2 and for d/R, \ P /R S> 1. 

(b) Interaction involving a Drude (D) cylinder with another 
Drude cylinder, a plasma cylinder, a perfect metal cylinder or 
a plate. The separating curve is given, up to logarithmic cor- 
rections, by d/R ~ \fd/\r. The shown expressions hold for 
d/R, d/\ a 3> 1 and d 2> Ap/A CT . Depending on the relative 
size of length scales, different regimes can be reached: Arrow 
(1) corresponds to an increasing distance d which ultimately 
leads to a non-universal interaction. Arrow (2) indicates an 
overall increase of the geometry (i.e., d/R fixed) with constant 
conductivity leading to a universal interaction. 



£' <C 1 or d <C R 2 /X a , we reproduce the universal (per- 
fect metal) asymptotic energy of Eq. ( 14 ). In the opposite 
limit £' 3> 1 or d 3> i? 2 /A CT , the asymptotic energy reads 
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16X a d 3 ln(2d/R) 



(20) 



Similarly, for a Drude cylinder with the characteristic 
length A CT parallel to a plasma (or perfect metal) cylinder, 
in the limit of d 3> A^ and d ^S> X 2 /X a the asymptotic 
Casimir energy is 

£ ir 2 hcR 2 



\n(X a d/R 2 ). 



(21) 



L 16X a d 3 ln(2d/R) 

These two limiting cases for £' are related to two differ- 



ent scaling regimes that are separat ed, up to logarithmic 
corrections, by the curve d/R ~ ^d/X a , see Fig. |l[b). 
The unconventional feature corresponds to the fact that 
the interaction is universal at shorter distances where 
d <C R 2 /X a . If the distance is increased beyond this 
crossover scale (with all other length scales kept fixed, see 
arrow (1) in Fig. [lib), the interaction becomes material 
dependent and, up to logarithmic corrections, scales as 
R 2 /(X a d 3 ) for a Drude cylinder interacting with another 
Drude or a plasma or a perfect metal cylinder. How- 
ever, if the radii of the cylinders are increased in the 
same way as their distance (d/R fixed, see arrow (2) in 
Fig. Pub)), finite conductivity becomes unimportant at 
large distances and the interaction assumes the univer- 
sal form. An intuitive explanation of this non-universal 
large distance behavior is given below. It is important to 
note that all forms of these metallic interactions decay 
much slower than the Casimir energy of two insulating 
cylinders which for d ^$> R scales as hcLR 4 /d 6 with a 
material dependent coefficient. 



B. Cylinder parallel to a plate 

In this section we consider a cylinder with radius R 
parallel to a plate. We show, similar to parallel cylinders, 
the existence of two different scaling regimes that are 
separated by curves given by the same expressions that 
we found for two cylinders, see Fig. IT] In order to find the 
asymptotic large distance interations, we employ again 
Eq. ([9]). The trace of the matrix N in Eq. (6l) in the limit 
of large separation d 3> R is approximated by 



Tr[N^]=^ 



EEkkEE 
l 00 ' 



(22) 



Ijjf = 1 and 



Note that for perfect metal plates Tjf 
for small k at fixed kj_/n one has for the plasma model 
T k ± = T k{ = l+0(\ p n) and for the Drude model T%_ = 
T^ = 1 + 0(\„k). Therefore, at large distances the 
material desription of the plate is unimportant and to 
leading order in R/d one gets 



Ifcf 



dk v =K Q (2pd). (23) 
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Using Eq. p2| in Eq. Q, and changing again considering 
k = pcos(9)/d and k z = psin(8)/d, we obtain 
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n "irH^r < 24 > 



where the functions C(p,6) = Ci(p,6) are given by the 
expressions below Eq. |TT| ) 

For a perfect metal cylinder, the integrals can be cal- 
culated in a straight forward manner, resulting in the 
universal energy 
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For a plasma cylinder with the plasma wavelength 
A p , after performing the radial and angular integrals in 



Eq. ( 24 ) , we obtain 
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(26) 



with g(x) = (1 + x) 1 ' 2 . As in our analysis for paral- 
lel cyli nder s we consider two different cases for £ given 
in Eq. jTl. If £ < 1 or similarly \n(2d/R) > A^/i? 2 , 
we reproduce the universal energy, Eq. (25 1. In the op- 



posite limit of exponentially large distances ^ > 1 or 
d ^$> Rexp(X 2 /R 2 ) the Casimir energy is non-universal 
and we obtain 
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(27) 



For a Drude cylinder with the characteristic length X a 
parallel to a metallic plate, in the limit f « 1 or 
d <C R 2 /X a , the integrand in Eq. (24) becomes indepen- 



dent of 9 and we reproduce the universal Casimir energy 
in Eq. (25). Hence, similar to the case of two cylinders, 



the interaction approaches a universal form below a ge- 
ometry and material dependent crossover distance. This 
counterintuitive result shall be discussed below. In the 
opposite limit (' > 1 or d > R 2 /X ai the asymptotic 
energy becomes non-universal and reads 



£ _ irhcR 2 \n(X (7 d/R 2 ) 
L ~ 8X~7d? 



(28) 



Based on the studies described above, we conclude that 
the Casimir interaction between a metallic cylinder and 
a plate decays slower than that between an insulating 
cylinder and a plane for which the energy scales as £ ~ 
HcLR 2 /d 4 for d^> R [UJ- 



IV. NUMERICAL RESULTS 

In this section, we compute the Casimir energy based 
on Eq. ([5]) at zero temperature. Our results are obtained 
by numerical computation of the determinant and the 
integrals over k and k z . Note that for a cylinder parallel 
to a plate in addition to the k and k z integrations, one 
has to compute the integral over k y for each element of 
the matrix N, see Eqs. ([6]) and Q. The matrix N (and 
hence the sum over n in Eq. p])) are truncated at a finite 
partial wave number n ma x- 

We chose n max such that the result for the energy 
changes by less than 0.01% upon increasing n max by 10. 
The required value of n max diverges when the surface-to- 
surface separation h between the objects (for cylinders 
h = d — 2R and for a cylinder and a plate h = d — R) 
tends to zero. For example, for 0.6 < h/R < 1.0, we used 
n max = 21, whereas for h/R = 0.6 and 0.5, one needs 
"-max = 31. For h/R — 0.05 we set the value n max = 191. 
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FIG. 2. Ratio of numerically computed energies i? for re- 
alistic metals and -E pcr f for perfect metal against h/R. (a) 
two identical cylinders, (b) a cylinder and a plate. The pa- 
rameters for the plasma model are X p /R = 0.05, 0.5 (open 
circles and open squares respetively) . For the Drude model 
the same values of X p /R are used; X p /R = 0.05, 0.5 (filled and 
open triangles, respectively) and X p /X a = 27.4. Insets: Short 
distance range. 



To reach sufficient numerical accuracy in the computa- 
tion of det(l — N) we have computed the Bessel functions 
with quadruple precision and employed similarity trans- 
formations for N by using the DEGBAL routine of the 
LAP A CK library with quadruple precision [25J . 

In Fig. [2j we show our numerical results for two par- 
allel cylinders and also for a cylinder parallel to a plate. 
The graphs show the Casimir energies for the Drude and 
plasma cylinders, normalized to the energies of perfect 
metal cylinders. For the numerics we used X p /R = 0.05 
and 0.5 with X p /X a = 27.4, corresponding to gold for 
which A p = 137 nm and X a « 5 nm. Figure [2] clearly 
shows the material dependence of the Casimir energies. 
At large separations, for the plasma model, the ratios 
of E/E per f approach one. This is due to the fact that 
Xp/R < 1 and we are in the universal regime, see Fig. [I] 
For Drude cylinders, the quantity £' determines the be- 
havior of the curves. At X p /R = 0.5, one has 5 x 10~ 4 < 
£' < 0.06 and for X p /R = 0.05, 5 x 10" 5 < f < 0.006 for 
the range of h shown in Fig. (pi). Since £' <C 1, we ex- 



1.00 
0.95 
0.90 

i 0.85 

U3 

£}" 0.80 
0.75 
0.70 



• ■ 
• ■ 
• ■ 



° oS oSSS8 



• ■ 
• ■ 
• ■ 
• ■ 
• ■ 



ll 



0.65 L- 
0.0 



0.2 



0.4 0.6 

h/R 



0^8 



FIG. 3. Ratio of the numerical results for the Casimir energy 
shown in Fig. [2] and the PFA energy based on the Lifshitz 
theory for the plasma model with X p /R = 0.05 (squares) and 
X p /R — 0.5 (circles). The ratio is shown as a function of 
the surface-to-surface distance h. Filled and empty shapes 
representing the data for parallel cylinders and cylinder-plate, 
respectively. 



pect from our asymptotic computations that the Casimir 
energy is close to the energy for perfect metal cylinders. 
Fig. [2] indeed shows a plateau at intermediate distances 
that is approaching the perfect metal energy E per f. This 
approach is better for X p /R = 0.05 which corresponds to 
a smaller £'. At small distances, none of our asymptotic 
results applies and the actual energy is more strongly re- 
druced compared to E per f . With increasing distance, we 
expect at £' ~ 1 a crossover to the non-universal asymp- 
totic energy of Eq. ( 20 1 . While this crossover is not fully 
shown in Fig.[2j the descrease of the energy ratio E/E per f 
with increasing separation is a precursor of this crossover. 
The same arguments apply to the interaction of a Drude 
cylinder with a plate. 

We now compare our numerics with the PFA results at 
short separations and with the asymptotic results at large 
separations for both the plasma and the Drude models. 



PFA versus numerics 



The PFA energy is obtained by integrating the PFA 
force F = 2tt RE p \ &tes {h) with respect to h, where 
E p ia,tes{h) is the energy of two parallel plates at distance 
h given by the Lifshitz formula [29 using the dielectric 
function of Eq. pi). Fig. [3] shows the numerically com- 
puted Casimir energy for X p /R = 0.5 and 0.05, normal- 
ized to the PFA energy for parallel cylinders and a cylin- 
der parallel to a plate. We find similar results for the 
Drude model. The energies associated with the Drude 
model are not shown here since they collapse on the data 
for the plasma model at short separations. Our data sup- 
port the consistency of the PFA in the limit of vanishing 
separations. 
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FIG. 4. Two identical cylinders: Ratio of the numerically 
computed energy for the plasma model (a) and the Drude 
model (b) and the corresponding universal (blue solid line) 
and non-universal (red dashed line) asymptotic results with 
Xp/R = 5 and X a = A p /27.4. The universal asymptotic en- 
ergy is g iven by Eq. ( 14 1 and non-universal ones are given 
by Eqs. (17 1 and (20 1. The yellow strip shows the crossover 



FIG. 5. Same as Fig. [4] but for a cylinder and a plate with 
the universal asymptotic energy given by Eq. ( 25 1 and non- 
universal ones by Eqs. (1271) and J28|. 



region between the Drude and plasma asymptotic energies. 



B. Asymptotics versus numerics 



Figures [4] and [5] show the ratio of the computed ener- 
gies and the corresponding asymptotic results (universal 
and non-universal regimes) versus ln(2d/R) for two cylin- 
ders and a cylinder parallel to a plate, respectively. The 
parameter that couples shape (radius) and material prop- 
erties is chosen as X p /R = 5. Figures Sa) and[5|a) show 
that for plasma cylinders, at intermediate separations, 
the energy normalized to the non-universal asymptotic 
energy is approaching unity whereas at asymptotically 
large separations the energy normalized to the univer- 
sal asymptotic energy is tending to unity. On the other 
hand, Figs. Kb) and p>tb) show that for Drude cylin- 
ders, at intermediate separations, the energy normalized 
to the universal energy is approaching unity whereas at 
asymptotically large separations the energy normalized 
to the non-universal energy is tending to unity. These 
figures confirm the validity of the crossover regime shown 
inFig.[T} 



V. SUMMARY 

In summary, we have calculated the Casimir force be- 
tween two metallic cylinders and a metallic cylinder par- 
allel to a plate. The energy is calculated numerically for 
a large range of separations. We also find asymptotic en- 
ergies for large separations and confirmed their validity 
with the numerical results. Furthermore, we showed that 
the numerics tend to the PFA energies at short separa- 
tions. 

The interesting phenomenon in our results is that the 
Casimir interaction involving Drude cylinders approaches 
a universal form of interaction at intermediate separa- 
tions and becomes non-universal (material dependent) 
at larger distances d > R 2 /X a . This behavior can be 
explained in terms of the size of the collective charge 
fluctuations in a Drude metal. However if A p <C R then 
plasma oscillations are supported by the cylinders and at 
asymptotically large separations the interaction energy 
docs show universality. 

Based on recent experiments, the interactions between 
metals might not be consistent with the Drude model 
|30) . The asymptotic energies that we found in this work 
can be used to provide a clearer distinction between the 
Drude and plasma model predictions as compared to two 
plates or a plate and sphere [16]. An estimate of the 
interaction between two gold cylinders with R = lOnm, 
length L = lOO/im, X p = 137nm and A CT = 5nm at a 
distance d — 200nm yields a force of w lpN within the 



plasma description and ~ 27pN within the Drude model. 
These forces are in the experimentally detectable regime. 
The significant feature of the interaction between a 
Drude cylinder with another Drude cylinder or a plate 
is that upon increasing the separation, the interaction 
can move from a universal regime to a non-universal one. 
This behavior can be understood from the wave equa- 
tion for the electric field inside a Drude cylinder. For 
imaginary frequencies u> = ion, the Helmholtz operator 
V 2 + e(uj)(cj/c) 2 for a good Drude conductor becomes 
V 2 — 8Tr 2 K/X a . We are interested in the maximal wave 
length of the field and hence charge fluctuations for a 
given k. With the smallest transverse wave vector k x , 
k y ~ 2tt/R we find the dispersion relation 



|fc 2 |~irV«/Kc-i, Kc = \ a /R 2 



(29) 



Hence, collective charge fluctuations on arbitrarily large 
scales exist only for k > k c which is a consequence 
of dimensionality that does not appear in the absence 
of transverse constraints (R — > oo). For k < k c 
charge fluctuations break up into clusters of typical size 
~ R/y/l — k/k c due to finite conductivity. The spec- 
tral contribution to the interaction between cylinders at 
distance d is peaked around k ~ 1/d. If d < 1/k c 
(d/R < y/d/\ a , see Fig. |l[b)), collective charge fluctu- 
ations contribute strongly to the interaction and render 
it universal similar to perfect metal cylinders for which 
k c ~ 1/cr — > 0. In the asymptotic regime with d > l/n c 
(d/R > y/d/Xu, see Fig. Illb)), finite conductivity pre- 
vents fluctuations on arbitrarily large scales and hence 
the interaction is proportional to a, i.e., non-universal. 
It is important to note that as R goes to zero, k c be- 
comes larger, and in consequence the finite conductivity 
of the cylinder becomes more important. 
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Appendix A: T-matrices 
1. T-matrix of a cylinder 



the fields inside and outside the cylinder by matching 
the boundary conditions at the cylinder surface. Finally 
in part (Ale) we show that the derived T-matrix in the 
limit of perfect conductivity (e — > oo) agrees with the 
T-matrix of a perfectly conducting cylinder. 



a. Vector cylindrical harmonics 

Vector cylindrical harmonics provide a basis in which 
divergence-less solutions of the vector Helmholtz equa- 
tion can be expanded. As in the spherical case, the vector 
harmonics can be obtained by applying the curl operator 
to a vector field $fc z „(ft,x) that is given in terms of the 
scalar harmonics. In the cylindrical case, however, this 
construction is simpler than in the spherical case since 
the vector field $fc a „(«, x) can be chosen to be parallel 
to the z-axis and hence does not change direction as it 
does in the spherical case where $; m («, x) ~ x. Depend- 
ing on the polarization we defined the two sets of regular 
cylindrical harmonics 



M£ B (*,x) 

N£(«,x) 



1 



\fk~i - ;•■' 



2 _|_ K .2 
1 



K \A1 



K 



Vx$£ g >,x), (Al) 

VxVx<^ s >,x),(A2) 



for magnetic multipoles (TE waves) and electric multi- 
poles (TM waves) , respectively, with p = \/k 2 + k 2 and 
the vector field 



*&(".*) 



zl n (pp)e 



/ nV ik z z 



(A3) 



where p, 9 and z are the cylindrical coordinates of x. The 

analog basis M£ u ^(k,x), N£ u ^(k,x) for outgoing waves 

is obtained by replacing /„ by K, : 

waves, i.e., 

lations M^ 

explicit form, they read 



VML„ = VN| a „ = 



These are transverse 
. They obey the re- 



^V x M? 



In 



MS>,x) 



VP 



'-In(PP)p-l'n(PP)Q 



e ine e ik ' z , (A4) 



N r fe ? n («,x) 



ik z l' n (pp)p 
-pl n (pp)z 



ik. 



In{p P )0 



VP 



It is analogous for the outgoing waves. 



In this subsection, we derive the T-matrix of a dielec- 
tric cylinder that is placed in the vacuum. For this pur- 
pose, in part ( A 1 a I we find a solution to the vector wave 
equation in terms of vector cylindrical harmonics. In 
part (Alb) we expand the electromagnetic field inside 



and outside of a cylinder in the basis of the solutions 
presented in the previous part. Then we find the expan- 
sion coefficients, which are the T-matrix elements, for 



b. Scattering amplitudes 

We consider an infinitely long dielectric cylinder with 
e(icK), p(icK) and radius R in vacuum. We expand the 
electromagnetic field inside and outside the cylinder in 
the bases of Eqs. ([Al]), ( |A2[ ) and the corresponding bases 
for outgoing waves. The expansion coefficients for the 
field inside and outside (T-matrix elements) follow from 



the matching conditions at the cylinder surface for the 
field components that are parallel to the surface. 

For an incident magnetic multipole (TE) field, we make 
the scattering ansatz for the electric field modes 



E£„(k, x) = MlJ n (n, x) + T^fMTXiK, x) 

+ T£fN££(*,x), (A6) 

outside the cylinder and 



E 



fc , n («,x) =<-M^„( K ,x) + A^N^ n («,x) ) (A7) 



inside the cylinder where M]?^(«,x), NT^(«,x) are 
given by Eqs.. (Al), (A2) with k replaced by 



ft a/ e(icK,) n(icn) . For an incident electric multipole (TM) 
field, the ansatz becomes 

E£„(«,x) = N^(M)+C M S(M) 

+ T^Nr t re ( K ,x), (A8) 



outside the cylinder and 






,(*) 



Af^M r ^ m («,x) 



Af^N^^x), (A9) 



inside the cylinder. 

The continuity conditions require that the tangen- 
tial fields E z , Ej,, H z — B z /fi and H^ = B^j\x are 
continuous across the cylinder surface. Using the 
explicit expressions of Eqs. ( |A4 ), ( |A 1 a ) these conditions 
lead for each type of multipo e fields to a set of four 
linear equations for the expansion coefficients. Using 
B = — (1/k)V x E and setting p' = ^/e/iK 2 + k 2 these 
equations can be written for incident magnetic (TE) 
waves as 



M, 



k z n 



( A h™\ 

rpMM 

k*n 

aME 

^k,n 



( E MpR) \ 

I'JpR) 

o 

\$tMpR)J 



(A10) 



with the matrix 



Mi 



k z n 



( fjniP'R) 

KXp'R) 
o 



lK n { P R) 
-K( P R) 







'" -MP'R) -%K n {pR) 



JtMp'r) 



- E K K n (pR) 



Kwk-MP'R) -^K n ( P R) y/e/nI> n (p>R) -K(pR) ) 



(All) 



For incident electric (TM) waves the linear equations and 



are 
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( A i M A 

rpEM 

k z n 
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( k ° \ 

IUpR) 
V I'nipR) J 



(A12) 



with the same matrix Mk z n as before. The solution to 
these equations for the T-matrix elements can be ex- 
pressed as 



Ai = 
A 2 = 

A 3 = 

A 4 = 



I' n (p'R) 1 K' m {pR) 



p'RI n (p'R) e P RK m (pR) 

r n (p'R) 1 K' n (pR) 

p'RI n {p'R) npRK n {pR) 

r n ( P 'R) i r n {pR) 

p'RI n (p'R) epRInipR)' 

rjj/R) i rjpR) 

p'RI n {p'R) i-ipRI n {pR) ' 



(A17) 
(A18) 
(A19) 
(A20) 



rpME 
-Lhn 



pAIM 
k z n 



rpEE 

k z n 



rpEM 

k z n 



I n (pR) A 1 A 4 + K 2 
K n (pR)A 1 A 2 + K^ [ ° j 

I n ( P R) A 2 A 3 +K 2 
K n {pR)A 1 A 2 +K 2 ' l } 

K 1 

^IJl(pR) 2 K n (pR) 2 A 1 A 2 + K 2 ' 

(A15) 



with 



Notice that in general the polarization is not conserved 
under scattering, i.e., Tf * f ^ ^ Tjf%. 



c. Limit of perfect conductivity 

We consider the limit e —¥ oo with (i fixed. Then p' — > 
ejjK and K ~ 1/ ' JT\i. In addition, we have for large p' 



K 



nk 7 . 



e^R 2 K \p' 2 p 2 
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I'nJP'R) 1 

p>RI n (p>R) pfR 
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so that 



Ai, A 3 

A 2 



1 



e/iRk 

1 K' n [pR) 
(ipRK n (pR) 



+ 0{e- 1 ), (A22) 

+ 0(e- 1 / 2 ), (A23) 

(A24) 



HpRI n {pR) 



This asymptotic forms show that the T-matrix elements 
that couple TM and TE waves vanish as 



T, 



ME 
k z n 



1 



0. 



(A25) 



in the limit e — > oo. The T-matrix elements that couple 
like polarizations simplify substantially. Since for e — > oo 



A!A 4 -jy 2 ; I' n {pR)K n (pR) 
A 1 A 2 -K* I n (pR)K' n (pR) 
A 2 A 3 - X 2 



1 



AiA 2 - K 2 
we get the simplified expressions 

I'nipR) 



T, 



MM 



T, 



K' n {pR)' 

EE _ I n (pR) 



k z n 



K n {pR) ' 



rpEM _ rpME _ Q 

k z n k z n ^* 



(A26) 

(A27) 

(A28) 

(A29) 
(A30) 



It is easily checked that these are the T-matrix elements 
for a scalar field with Neumann boundary conditions 
(magnetic or TE modes) and with Dirichlet boundary 
conditions (electric or TM modes). Hence, in the limit 
of perfect conductivity the EM scattering problem for a 
cylinder separates into two independent scalar problems, 
one with Dirichlet and one with Neumann boundary con- 
ditions. 



2. T-matrix of a plate 

The T-matrix elements of a plane is given by its Frensel 
coefficients [TT] 



7£ =r-(z CK ,(l + ki/^)-^), 



(A31) 



with P the polarization index, kj^ the momentum per- 
pendicular to the x direction, and r® the Frensel coeffi- 
cients 

V ' p,(icK)+^l + (n 2 (icK,)-l)x 2 V ; 

r E (icK x) = e (' CK )~v /I +(" 2 (' CK )- 1 ) x2 ( A33 \ 

V ' ; £(icK)+ v /l + (n 2 (icK)-l)K 2 ' V ' 

here e, /x and n(icK) are the dielectric response function, 
magnetic permeability, and the refractive index of the 
plate, respectively. The refractive index in terms of the 



dielectric function and the magnetic permeability of the 
plate is given by 



n(icn) = \J e(icK) /j,(icK) . 



(A34) 



The T-matrix of a dielectric plate Eq. ( A31 ) is diagonal 
with respect to the polarization indices. 



Appendix B: Translation matrix 

According to Graf's addition theorem, the following 
relation holds, 

K m (p rj )e lm ^ = 

OO 



(Bl) 



where r, 



Rji with r,: 



i), Rj 



Rji(cos(f>ji, sin 4>ji) are two-dimensional vectors in the 
xy-plane. We consider translations in 3D that are per- 
pendicular to the z-axis with the translation vector X^^ = 



(Rji,0), i.e., we set Xj 



, Zi) with Zi 



Zj. Since the 



curl operator commutes with translations, from the def- 



initions in Eqs. (Al), (A2| and the addition theorem of 



Eq. (Bl) follow the translation formulas from outgoing 



to regular waves 



M° fc f m ( Xj -)= £ ^„ m M^( Xi) (B2) 



n— — oo 



Kl?m(Xi)= E U ZnmKU*i) (B3) 



(B4) 



n— — oo 

with the translation matrix 

Kn m =K m -n(p\X Jt \)e-^~ n ^. 

From this we make two important observations: Trans- 
lations conserve the polarization, i.e., they do not couple 
magnetic and electric modes, and the translation matri- 
ces are diagonal in k z . The conservation of polarization 
leads to a diagonal translation matrix that acts on the 
full set of electric and magnetic modes. 



Appendix C: Conversion matrix from vector plane 
wave basis to cylindrical vector wave basis 

In this section we show the conversion matrix elements 
from Ref . [TT] . The cylindrical vector wave functions are 
given in Eq. (Al) and (A2), which decay along the — z. 



We consider regular vector plane wave functions that de- 
cay along the — x axis 
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yj 'K, 2 +k^+k^x+ik v y+ik^Z'- 



x. (C2) 



The vector plane wave functions can be written in terms 
of vector cylindrical wave functions, 

M rc V( K ,x) =J2 D k»nM,k ± MM Te Z kzn ( K ,x) 

n 

+D kznEM±M W C h z n(K,x) , (C3) 

N™\ ± (K, X) - J2 Dk,nM,k ± EM Te h*n(K, x) 



+-Dfc 2 nE,kx£;N r0g fez „(K, x) , 



(C4) 



using the conversion matrix elements 
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k z nM,k±M ~ ~ l 



fc2 + k 2 



D kznE ,k ± M = i , £, (Vi + e 2 + e 



ki + ki 



Dk z nE,k ± E — Dk x nM,k±M i 
Dk z nMM±E = —Dk z nE,k±M, 



where £ 



y^+fef 



and k ± = (k y ,k z ). 
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